#pragma GCC optimize(2)
#include <iostream>
#include <queue>
#include <map>
#include <unordered_map>
#include <vector>
#include <algorithm>
#include <cmath>
#include <string>
#include <cstring>
#include <set>
#include <unordered_set>
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
const int mod=9973;

void mul(ll* a,ll* b,const int row,const int col)
{
    ll c[row][col];
    memset(c, 0, sizeof(c));
    for(int i=0; i<row; i++)
        for(int j=0; j<col; j++)
            for(int k=0; k<col; k++)
                c[i][j] = (c[i][j] + a[i*col+k]*b[k*col+j] % mod) % mod;
    memcpy(a,c,sizeof(c));
}

void qpow(ll* a,ll* b,int n,int row,int col)
{
    while(n)
    {
        if((n&1)==1)
        {
            mul(a,b,row,col);
        }
        n>>=1;
        mul(b,b,row,col);
    }
}

int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n,p;
        cin>>n>>p;
        ll a[n][n];
        ll b[n][n];
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                cin>>b[i][j];
                a[i][j]=0;
            }
        }
        for(int i=0;i<n;i++)
        {
            a[i][i]=1;
        }
        ll* pa=&a[0][0];
        ll* pb=&b[0][0];
        qpow(pa,pb,p,n,n);
        int ans=0;
        for(int i=0;i<n;i++)
        {
            ans=(ans+a[i][i]%mod)%mod;
        }
        cout<<ans<<endl;
    }

    return 0;
}
